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Thermal-Capillary Model 


• TDMUTOM FIELD DETERMINED 
CONDUCTION DOMINATED HEAT TRANSFER 
IN ICLT AND CRYSTAL. 

• WLT/SOLID INTERFACE SHAPE DETERMINED AS 

MELTING POINT ISOTHERM. 

• ICNI SOUS SHAPE DETERMINED BY BALANCE OF SURFACE 
TENSION AND HYDROSTATIC FORCES. 

• THICKNESS OF SHEET DETERMINED BY CONDITION 
FOR EQUILIBRIUM GROWTH ANGLE. 

• DISTANCE FROM DIE TIP TO HEIGHT OF HF.LT POOL 
NETS RIPEN KNCI PRHVIIMN IN iKNISCU* REGION. 
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SILICON SHEET 


Model Description 

• Find* Element Analysis to Provide: 

- Coupled Solutions for Heat Transfer and Capillarity 
In Three-Phase Domain of EFG Die/Melt/Crystal. 

- Three Unknown Boundaries for Crystal Thickness, 
Melt/Solid, Melt/Gas Interfaces. 

• Two-Dimensional Navler-Stokes Flow Field 

in Die Top and Meniscus (Interface) Melts 

• Diffusion Equation Solutions 

for Segregated Aluminum Dopant 


Physics Is Best Explained for Growth Into Ambient 
at Uniform Surrounding Temperature 
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CRYSTAL THICKNESS (cm) 


ORIGINAL PAGE fS 
OF POOR QUALITY 


SILICON SHEET 


Operating Diagram for EFG 



GROWTH RATE (cm/mi n) 


INLET MELT TEMPERATURE, t d 

• MELT BECOMES SUPERCOOLED FOR 
LOU DIE TEMPERATURES 


Crystal Thickness Is Maintained Constant 
by Simultaneously Changing v and Tp 



T D CK) 

V(cm/min) 

"Vc, 
"3y” cm 

1773 

0.96 

875 

1763 

1.81 

743 

1753 

2.57 

612 

1743 

3.41 

475 

17.' 

4.28 

333 

1723 

5.15 

186 

1713 

6.11 

31 


AXIAL DISTANCE (CM) 


• A MAXIMUM GROWTH RATE EXISTS WHERE TEMPERATURE GRADIENT AT INTERFACE DROPS TO ZERO. 
ONSET (*• SUPERCOOLING. 
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SILICON SHEET 


Ambient Temperature Distribution Used 
for Comparison With Experiments 



Comparison of Calculations and Measurements for 
Prediction of Thickness Variation With Pull Rate 



MELT-SICE INTERFACIAL 
TEMPERATURE GRADIENT (k/cm) 






ORIGINAL PAGE « SILICON SHEET 

OF POOR QUALITY 

Operating Region Predicted by Finite Element Analysis 



Dopant Segregation Effects 

Distribution ot Dopant Through Ribbon Thicknass Depend* on: 

(1) Interface Shape - Heat Transfer Solution 

(2) Solidification Flow Field - Navier Stokes Solutions 
Parameterized by Die Geometry. Meniscus Configuration 

(3) Surface Tension-Driven Flow 


Accurate Calculation of Dopant Segregation 


• r in ITS ELEMENT ANALYSIS Or THERMAL-CAPILLAR! MODEL SETS 
SHAPE or WLT. 

t VELOCITY riELD IE *LT CO# UTED SX PINITE ELOCVT 
SOLUTION rULL NAVIE1 -STOKES EQUATIONS IN «LT. 

• DOPANT CONCENTRATION FIELD CALCULATED BT 
SOLVING SPECIES CONSENVATION EQUATION 

V* e - Pe* V- (o y) • 0 

- 213V* / 0 


• NO ADJUSTABLE PAAAWTERS IN CALCULATION I 



579 



SILICON SHEET 


Lateral Dopant Segregation 



LATERAL DOPANT SP ACRE GAT ION PREDICTED BY CALCULATION OF DETAILED VELOCITY 
AND CONCENTRATION FIELDS IN THE KELT USING HEAT TRANSFER ANALYSIS. 
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SILICON SHEET 


Aluminum Distribution Across Thickness of Si Sheet 
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Surface-Tension-Driven Flow 


v 


_ do AT* 

"df 7 


For alHcow ■httl Sf Q; 
AT* ~ 5-20°K (modal) 


do 

dT 

M 

V 


>0.2 dynas/cm-K (Hardy) 

0.0080 polaa 

120 cm/aac » V, (0.04 cm/aac) 




SILICON SHEET 


Summary 


• Fit of Finito Element Model Prodictions to Experimental Data Achieved For 

- t-V* Relationship 

- Growth Rate Limits 

- Aluminum Dopant Segregation Dependence on t, V # 

With One Adjustable Parameter 

• Comprehensive Nature of Finite Element Model Demonstrated 

for Predicting Process Variable Relationships, 

Interface Configuration and Dopant Segregation 


NONLINEAR ANALYSIS OF MORPHOLOGICAL TRANSITIONS DURING 
DIRECTIONAL SOLIDIFICATION TRACES EVOLUTION OF PLANAR 
INTERFACE INTO HIGHLY STRUCTURED CELLULAR AND 
DENDRITIC FORMS. 



Plans ceils dendrites 


• Analysis combined finite-element-methoos for solving 

EQUATION-SET FOR TRANSPORT AND MORPHOLOGY WITH NEW 
METHODS IN NONLINEAR ANALYSIS 


e First to predict transition to deep cells 


• SHOW FUNDAMENTAL MODE FOR REDUCTION OF CELL WAVELENGTH 

e First to explain role of grain boundaries in onset op 

CELLULAR STRUCTURE 

• Analysis capable of prediction of onset of dendritic 
growth 
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SILICON SHEET 


Two-Dimensional Model of Interface Morphology 


Malt 
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Models for Studying Morphological Structure 


MODEL 

LATENT 

HEAT 

THERMAL 

CONDUCTIVITIES 

CONVECTIVE 
heat transport 

SOLIO 

DIFFUSION 

REFERENCE 

THERMAL-SOLUTAL 

(T3M) 

YES 

NOT EQUAL 

YES 

YES 

UNGAR ET Al_ U9B4) 

ONE-SIDED TSM 

YES 

NOT EQUAL 

YES 

NO 

MULLINS ANO SEKERKA 
4 Appl. PNys. 34 323(1989) 

EQUAL 

CONDUCTIVITY 

YES 

EQUAL 

YES 

NO 

MC FADDEN ANO COR 1 ELL 
PKytico D in pr«» (1984) 

SOLUTAL MODEL 
(SM) 

NO 

EQUAL 

NO 

YES 

UNGAR AND BROWN, 
Phyr Rn.B. (1904) 

symmetric sm 

NO 

EQUAL 

NO 

YES 

l 4* 

LANGER. Rev. Mod. Phft. 
52, 1 (I98T5T 

ONE-SIOED SM 

NO 

EQUAL 

NO 

NO 

UNGAR ANO BROWN 
Phy*. Rt*. B. 29. ISB7 (I9B4) 
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SILICON SHEET 



Progress Report 

\. FINITE ELEMENT PROGRAM FOR IDEALIZED GROWTH INCLINED GROWTH SYSTEM 
HAS JUST SEEN COMPLETED. RADIATIVE COUPLING BETWEEN MELT AND CRYSTAL 
SURFACES IN UNDERWAY. 

2. PROGRAM HILL BE READY FOR COMPARISON WITH EXPERIMENTS BY 01/01/1985. 

3. COUPLING OF THERMAL STRESS ANALYSIS TO THERMAL-CAP ILLARY NOOEL 
HILL BEGIN THIS FALL AMD BE COMPLETED IN LATE SPRING. 

A. FINITE CLEMENT ANALYSIS OF MICROSCOPIC INTERFACE MORPHOLOGY EETOND 
THE POINT OF LINEAR INSTABILITY HA if BEEN DEVELOPED. CALCULATIONS 
FDR SILICON UNDER SHEET GROWTH CONDITIONS WILL BE COMPLETED BY 
EAR! 7 S D RINC. 
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SILICON SHEET 


Sum nary 


Finite element analysis is being used on two length scales 

TO UNDERSTAND CRYSTAL GROWTH OF THIN SILICON SHEETS. 

1. Thermal-capillary models of entire ribbon- 
growth systems. Demonstrated for EFGj 
model presented for inclined-meniscus system. 

2. fllCROSCOPIC modeling of morphological structure 
of melt/solid interfaces beyond the point of 

LINEAR INSTABILITY. The FORMATION OF DEEP CELLS 
AND DENDRITES. APPLICATION TO SILICON SYSTEM IS 
UNDERWAY. 
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